Self-organized criticality in the Bean state in YBa 2 Cu 3 07_ ;E thin films. 
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The penetration of magnetic flux into a thin film of YBa2Cu307_ l; is studied when the external 
field is ramped slowly. In this case the flux penetrates in bursts or avalanches. The size of these 
avalanches is distributed according to a power law with an exponent of r = 1.29(2). The additional 
observation of finite-size scaling of the avalanche distributions, with an avalanche dimension D = 
1.89(3), gives strong indications towards self-organized criticality in this system. Furthermore we 
determine exponents governing the roughening dynamics of the flux surface using some universal 
scaling relations. These exponents are compared to those obtained from a standard roughening 
analysis. 
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The critical state in a type-II superconductor shows a 
powerful analogy to a granular pile, which was already 
noted by de Gennes in the 1960s [1]. With the advent of 
self-organized criticality (SOC) [2], the avalanche behav- 
ior of granular piles was intensely studied [3] . This is be- 
cause SOC was thought of as a general mechanism to ex- 
plain the intermittent behavior of slowly driven systems 
far from thermodynamic equilibrium. The experimental 
verification of SOC, however, was not straightforward. 
Power-law distribution of avalanches was not observed 
in many experiments [4] . One of the hall-marks of criti- 
cal behavior, finite-size scaling, was only found in a few 
cases, most notably in experiments on a one dimensional 
(Id) pile of rice [5], as well as in a Id pile of steel balls 
with a random distribution of balls in the bottom layer 
[6]. The only study of finite size scaling in a 2d system 
to date, is to our knowledge a study by some of us on the 
properties of a 2d pile of rice [7], which is qualitatively 
similar to the one presented here. 

Given the classical analogy with granular piles [8] , the 
critical state in superconductors was also quickly pro- 
posed as a SOC system [9], Experimentally, magnetic 
vortices are well suited to study SOC, since kinetic ef- 
fects, which can lead to deviations from critical behavior 
in sand-piles [5,10], are naturally suppressed due to their 
overdamped dynamics [11]. Just as in the case of granu- 
lar piles, however, the experimental confirmation of this 
conjecture has been controversial. While power-law be- 
havior in the avalanche distribution has been observed by 
most authors [12,13], finite-size scaling was not observed 
so far in the case of the critical state in superconductors. 
This is because most of the studies [12,14], were carried 
out using magnetization measurements, which only give 
information on the overall behavior of the whole of the 
sample. This corresponds to only considering off-edge 
avalanches in a sand-pile, which do not capture the full 
dynamics and may therefore give a flawed picture [10]. 
More recent investigations using arrays of miniature Hall- 
probes [13], do give insights into internal avalanches and 
the complete dynamics, but only give information from 



a few selected points in the sample, which makes testing 
for finite-size scaling impossible. 

In this letter, we study the local changes in the mag- 
netic flux over the whole central area of a sample. This 
is done via a highly sensitive magneto-optic setup, which 
can resolve flux densities of 0.2 mT over an area of ~ 
5x5 fim 2 . This implies that flux changes corresponding 
to 2.5 $o can be resolved, where $ = h/2e, is the mag- 
netic flux quantum (the flux of a single vortex). Fur- 
thermore, since we observe a large field of view, different 
size subsets can be studied leading to the possibility of 
observing finite-size scaling. Given the finite-size scal- 
ing exponents, one can determine the roughening and 
dynamical exponents of the flux-surface from universal 
scaling relations. Comparing these results with those de- 
termined directly from the surface properties [15], gives 
good agreement. Furthermore, the exponents fulfill the 
scaling relation a + z = 2, valid for the well known driven 
interface model given by the Kardar-Parisi-Zhang (KPZ) 
equation [16] (see below). 

The experiments were carried out on a thin film of 
YBa2Cu307_ x (YBCO), grown on a NdGaOs substrate 
to a thickness of 80 nm using pulsed laser ablation [17]. 
The pinning sites in the sample are uniformly distributed 
and consist mostly of screw dislocations acting as point 
pins [18] . A polarization microscope was placed in an Ox- 
ford Instruments cryomagnet capable of a maximum field 
of 1 T and cooling to a temperature of 1.8 K. The sample 
was cooled in zero field to 4.2 K, at which point the exter- 
nal field was slowly increased in steps of [1qA.H = 50 /iT, 
after which the sample was allowed to relax for 10 s before 
an image was taken. In a magneto-optical experiment, 
the flux density B z (x,y) at the surface of the sample is 
measured via an indicator layer showing a strong Faraday 
effect. The polarization of the incoming light is turned 
in proportion to the flux density in the indicator layer. 
A cross-polarized analyzer will thus admit light from re- 
gions with non-zero magnetic flux density [19]. However, 
in order to determine the rotation angle (and thus the 
flux density B 2 ) directly, including its sign, we employ 



1 



a recently developed magneto-optic image lock-in ampli- 
fier [20] , using a modulation of the incoming polarization 
vector. The output of the instrument gives directly the 
Faraday angle for each pixel, independent of spatial in- 
homogeneities in the illumination. 
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FIG. 1. Snapshot of a flux- avalanche. Shown is the dif- 
ference between two consecutive field- profiles AB(x,y). The 
scale goes from (white) to 0.4 mT (black). The area-integral 
over this difference corresponds to the size of the flux-jump 
A$. In order to determine the roughness exponent of the 
whole surface, we also need to determine the fractal dimen- 
sion of the active sites in such an avalanche (see text). The 
inset shows the flux distribution B z (x,y) in the same part of 
the sample. Here the scale varies from OmT (black) to 12mT 
(white). Lines of constant field can be identified by the dif- 
ferent shades. 

The data analyzed here come from a series of nine ex- 
perimental runs, each consisting of ~ 300 time-steps. Of 
these steps only the last 140 in each run were used, in or- 
der to have a critical state established in the whole region 
of the sample used for analysis. The size and shape of the 
avalanches was determined from the difference AJ5 z (x,y) 
of two consecutive images (see Fig. 1). From this differ- 
ence, the average increase in the applied magnetic field, 
due to the step-wise field-sweep, was subtracted in order 
to solely study the avalanches. The external fields were 
determined from a region well outside the sample, which 
was also in the field of view. Once the incremental field 
difference is determined, the size of an avalanche, corre- 
sponding to the displaced amount of flux A$, is obtained 
from AB via integration over the whole area 

s = A$ = i J \AB(x,y)\dxdy. (1) 

The average increase in the external field is in good agree- 
ment with the increase in flux density of 50 /iT between 
images. The resulting time series of the avalanche be- 
havior of all the experiments is shown in Fig. 2. In this 
figure, AB was integrated over the whole area of 180x180 



jiva 2 in each time step. As can be seen, the evolution of 
the magnetic flux inside the sample is intermittent with 
occasional large jumps. 
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FIG. 2. The time evolution of the magnetic flux inside the 
sample over all nine experiments. The magnetic field dif- 
ference AB has been integrated over the whole surface of 
180x180 fim 2 and the average flux increase has been sub- 
tracted. The evolution takes place in the form of flux jumps 
or avalanches of various sizes, which are summarized in the 
histogram of Fig. 3. The different experimental runs are sep- 
arated by arrows. 

In order to check the data for finite size scaling, we also 
integrated AB over subsets of the image of a linear size 
of L = 90, 45 and 15 yum respectively. The histograms 
of the avalanche size distribution for these data-sets are 
shown unsealed in Fig. 3a. As can be seen, the smallest 
avalanches correspond to a flux change of 2-3 &o, corre- 
sponding to the resolution of the measurement. Taking 
all of the data together we observe a power-law distri- 
bution over more than three decades. The slope of the 
black line gives the exponent of the distribution, r = 
1.30(5). In Fig. 3b we show the same data, but now the 
probabilities are scaled with s T and the avalanche sizes 
are scaled with L~ D . As can be seen, there is very good 
curve collapse indicating the presence of finite size scal- 
ing [21]. This means that the avalanche size distribution 
function is given by 

P(s,L) = s- T f(^), (2) 

where f(x) is constant up to a cutoff scale s co oc L D . The 
values of the exponents used to obtain curve collapse are 
t = 1.29(2) and D = 1.89(3). Note that here, we have 
carried out the finite size scaling by way of subdividing 
the whole image rather than carrying out experiments 
with different size samples. We have checked by means 
of simulations of the 2d Oslo-model that such a subdi- 
vision into finite-size samples leads to the same scaling 
exponents as a curve collapse of different simulations of 
finite size. Furthermore, an independent measurement of 
D, directly using a box counting method [22] in 3d yields 



2 



D = 1.92(5), 
scaling. 



consistent with the value from finite size 
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FIG. 3. (a) The direct avalanche size distribution for dif- 
ferent sizes of windows (180, 90, 45, and 15 /im respectively). 
There is an envelope to the distribution in the form of a 
power-law indicated by the black line over more than three 
decades with an exponent of r = 1.30(5). (b) The scaled size 
distributions showing a curve collapse. The exponents used 
in the finite-size scaling are r = 1.29(2) and D — 1.89(3). 



It has already been noted [23] that the front of pen- 
etrating flux shows kinetic roughening. Similarly, the 
surface B z (x,y) in two dimensions can be shown to be 
sclf-affine, with a roughness exponent a, characterizing 
the growth of the interface width w sa t with the system 
size, w sa t oc L a and a dynamic exponent z characterizing 
the saturation time of the width. In such an analysis of 
the present data, discussed elsewhere [15], the roughness 
exponent was found to be a = 0.73(5) and the dynamic 
exponent was found to be z = 1.38(10). Using the uni- 
versal scaling relations derived by Paczuski et al. [24] for 
various SOC models, these exponents characterizing the 
roughness and growth of the surface can also be deter- 
mined from the scaling exponents of the avalanche dis- 
tribution. This indicates the fact that the roughening 
of the surface originates from the avalanche distribution 



and its underlying dynamics. 

Let us first discuss the roughness exponent. Accord- 
ing to the finite size scaling, an avalanche of the cut-off 
size will scale like s co oc L D . Similarly, the size of such 
an avalanche will roughly be given by s co ~ w sa tL dB , 
where d,B is the fractal dimension of the area of an 
avalanche. Equating the two expressions for s co , one ob- 
tains w„ a t oc L D ~ dB and hence 



a = D — dg, 



(3) 



in agreement with Ref. [24]. In order to determine 
the roughness exponent from the avalanches, we there- 
fore have to measure the fractal dimension, cIb, of the 
avalanches area, which was done using a simple box- 
counting method [22]. To this end, avalanches which 
were one standard deviations bigger than average were 
binarized, yielding a distribution of active clusters used 
in the determination of the fractal dimension. The re- 
sult for one such cluster can be seen in Fig. 4, where 
the number of active pixels of the avalanche is shown as 
a function of the length scale. Averaged over all clus- 
ters analyzed, the fractal dimension, given by the slope 
of the line in Fig. 4, is ds — 1.18(5). From this we de- 
termine the roughness exponent as a = 0.71(5), in good 
agreement with that determined from a roughness analy- 
sis (via the correlation function and the power-spectrum) 
of the surface fluctuations [15]. 
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FIG. 4. Determination of the fractal dimension of the ac- 
tive clusters in an avalanche using the box counting method. 
In this determination, avalanches more than one standard de- 
viation larger than the average size have been studied. The 
analysis of one such cluster is shown. The fractal dimension 
is given by the slope of the straight line. Averaged over all 
clusters, this leads to a fractal dimension of the active clusters 
of d B = 1.18(5). 

The dynamic exponent, z, can be obtained from a sim- 
ilar argument. The scaling of the crossover-time with the 
system size is described by z, t x ex L z . This crossover- 
time can be roughly estimated from the time it takes for 
an avalanche of the cut-off size s co to appear. Since flux 
is added to the system in constant steps, the number 
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of vortices added until a cut-off avalanche occurs is given 
by <5<&i x - On the other hand, the flux added will also 
be moved in smaller avalanches. The total flux necessary 
to be introduced in order to obtain an avalanche of size 
s co can be estimated from integrating over all avalanche 
sizes up to the one of size s co : 

t x oc f sP(s)ds cx s\- T . (4) 
Jo 

Using s co oc L D , we obtain t x oc L D ( 2 ~ T ),which immedi- 
ately leads to the scaling relation 

z = D(2-t). (5) 

Again, this is also in agreement with the universal scal- 
ing relation derived by Paczuski et al. [24] . Inserting the 
values determined above, we obtain z = 1.34(4), again 
in good agreement with that determined via roughness 
analysis [15]. 

In conclusion, we have shown that the distribution of 
the size of flux jumps in an YBCO thin film is not only 
given by a power-law, but also shows finite size scaling. 
Due to universal scaling relations valid for a SOC state 
[24] , the exponents determined via finite-size scaling can 
also be used to describe other properties of the system. 
One such example is the roughening of the magnetic flux 
surface. Here the statistical properties of the self-affinc 
structure are built up by the penetrating flux, such as the 
roughness and dynamical exponent can be derived from 
the structure and dynamics of the flux-avalanches us- 
ing relations (3) and (5). These characteristic exponents 
can however also be determined directly, via a rough- 
ness analysis [15], the results of which can be compared 
to those obtained via the avalanche dynamics. As dis- 
cussed above, we find excellent agreement between the 
exponents determined in these separate ways. The crit- 
ical state observed in the YBCO thin films can be seen 
as a realization of a 2d roughening system, albeit with 
a self-organized dynamics. In this context, we note that 
the roughness and dynamic exponents fulfill the general 
KPZ scaling relation a + z = 2, which is an exact result 
also in 2d. 

In the future, the magnetic flux structures in super- 
conductors, as the ones studied here, may be used as 
an ideal experimental system with which to study non- 
equilibrium phenomena, especially those of granular mat- 
ter. In fact we find strong qualitative correspondence of 
the behavior of the vortices with that of a pile of rice 
[7] in terms of SOC behavior. However, in addition to 
the granular system, the superconducting system stud- 
ied here allows for the presence and control of quenched 
noise due to pinning sites, whose influence on the sta- 
tistical properties can then be studied experimentally as 
well. Moreover, the physics of the microscopic behavior 
of vortices is well studied [11], such that collective effects 
can be studied directly. 
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